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Abstract. A comparison between the oscillation frequencies of six multi-periodic 8 Scuti stars of the Pleiades 
cluster and the eigenfrequencies of rotating stellar models that match the corresponding stellar parameters has 
been carried out. The assumption that all the stars considered have some common parameters, such as metallicity, 
distance or age, is imposed as a constraint. As a result, we have a best fit solution associated with a cluster 
metallicity of [Fe/H]~ 0.067, an age between 70 x 10 6 and 100 x 10 6 yr and a distance modulus of mv — Mv = 5.60- 
5.70 mag. All the stars were found to oscillate mainly in non-radial, low degree, low order p modes. Estimates of 
mass and rotation rates for each star are also obtained. 
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1. Introduction 

The study of 5 Scuti stars is expected to provide an im- 
portant diagnostic tool for the structure and evolution of 
intermediate-mass stars at the main sequence evolution 
and thick-shell hydrogen burning stages. Having masses 
between 1.5 and 2.5 M@, these pulsating variables develop 
convective cores during their central hydrogen-burning 
phases that make them suitable for investigating the hy- 
drodynamic processes occurring in the core. 

Although most of the general properties of these stars 
are well understood, the application of seismological meth- 
ods to derive their properties has proved to be a diffi- 
cult task. The problems concerning the modelling of S 
Scuti stars have been discussed in details in recent re- 
views (e.g. Christensen-Dalsgaard 2000 Kieldsen 12001(1 . 
The most acute seem to be the mode selection mecha- 
nism and the influence of rotation. Concerning the for- 
mer, although there is fairly good agreement between the 
observed frequency range and that derived from instabil- 
ity calculations (e.g. Michel et al. I1999f) . far more modes 
than observed are predicted to be unstable. While it is ex- 
pected that the forthcoming asteroseismology space mis- 
sion COROT (Baglin et al. EHM)) will be able to disen- 
tangle the whole spectra of S Scuti stars, for the moment 
the mechanism that could cause an amplitude limitation 
in some modes is still unknown. As a consequence, mode 
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identification based on the comparison between theoretical 
and observed frequencies is difficult to attain in general. 
On the other hand, although several observational tech- 
niques for mode identification have been developed in re- 
cent years ( Watson IT5ggl Garrido et al. lT5501 Aerts 115551 
Kenelly & Walker H"55B1 Viskum et al. 11555)1 . in only a few 
cases have these techniques been successfully applied. 

It is well known that most S Scuti stars are rapid rota- 
tors. This has important effects on the mode frequencies, 
both directly as a consequence of the changes that must 
be introduced in the oscillation equations and indirectly 
through the changes in the equilibrium models. Therefore, 
successful analysis of rotating S Scuti stars requires that 
rotation should be taken into account consistently at all 
levels in the analysis. 

Despite these problems, in recent years several at- 
tempts have been made to interpret the observed spectra 
of S Scuti stars (Goupil et al. 119931 Perez Hernandez et 
al. 115551 Guzik & Bradley 115551 Pamyatnykh et al. 115581 
Bradley & Guzik 119961 Hernandez et al. 119981 Suarez et 
al. I2005[) . Here we address the problem of mode identi- 
fication for stars in an open cluster. These stars are ex- 
pected to have a common age and chemical composition. 
Furthermore, the distance to the cluster is usually known 
with high accuracy. The constraints imposed by the cluster 
parameters have proved to be very useful when modelling 
an ensemble of S Scuti stars. The best examples of such 
studies are found for the variables in the Praesepe cluster 
(Michel et al. 119991 Hernandez et al. [1998). In a similar 
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way, we consider here several S Scuti stars of the Pleiades 
cluster and search for a best fit solution in the sense of 
a set of stellar parameters that allows the simultaneous 
modelling of all the stars considered, and that satisfies all 
the observables, including the frequencies. 

The paper is organized as follows. In Sect. [21 we de- 
scribe the main characteristics of the Pleiades cluster and 
the six S Scuti stars under study. The modelling, the range 
of parameters and the calculations of the eigenfrequencies 
are discussed in Sect. El The details of the comparison be- 
tween observed and theoretical frequencies are discussed 
in Sect. The results and their discussion are given in 
Sect. [5] Finally, we present our conclusions in Sect.EJ 

2. The observational material 

The Pleiades (M45) is a young (~ 75-100 Myr), 
Population I cluster. Because of its proximity, the obser- 
vational parameters of the Pleiades have been intensively 
studied. In particular, the metallicity of the cluster is esti- 
mated to be between -0.14 < [Fe/H] < +0.13 (e.g. -0.034 
± 0.02 4 Boes gaard & FrielfTMOl 0.0 260 ± 0.103 Cayrel de 
Strobel lTWil -0.11 ± 0.03 Orenon l27)UTl) . 

Until recently, there was a dispute regarding the dis- 
tance of the Pleiades cluster from the Earth. While the 
determination of the cluster distance from direct par- 
allax measurements of HIPPARCOS gives 116.3±g | pc 
(my - My = 5.33 ± 0.06 mag, Mermilliod et al, ITM7|l . 
the previously accepted distance, which is based on com- 
paring the cluster's main sequence with that of nearby 
stars, was «130 pc corresponding to a distance modulus 
of about 5.60 mag (e.g. 5.65 ± 0.38 mag [Eggen lT^] . 5.60 
± 0.05 mag [Pinsonneault et al. ITM8] . 5.61 ± 0.03 mag 
[Stello & NissenEOnil). As will be discussed in Sect. [U 
the distance derived in the research presented here agrees 
with the latter value. 

To date, six S Scuti stars are known in the Pleiades. In 
a survey of variability in the cluster, Breger l|1972f) found 
S Scuti type oscillations in four stars. The remaining two 
were recently reported to be 5 Scuti stars by Koen l|1999[) 
and Li et al. ( 2002 ) . Some observational properties of these 
stars are given in Table ^ The projected rotational ve- 
locities, usini, were obtained from Morse et al. I|1991[) 
and Uesugi & Fukuda (1982). A 10% uncertainty in these 
quantities has been assumed. The spectral types (ST) and 
(3 parameters were obtained from the SIMBAD database 
in Strasbourg (France). The errors in the photometric pa- 
rameters (B — V) and my are estimated from the disper- 
sion between different measurements of these quantities. 
The errors in reddening are taken from Breger |I986J. We 
shall use these uncertainties in Section 3.2. 

In recent years, the 6 Scuti stars of the Pleiades cluster 
have been observed in several campaigns of the STEPHI 
multi-site network (Michel et al. 12000)). The information 
on the oscillation frequencies of the stars used here has 
been obtained from those campaigns and are published 
elsewhere (Belmonte & Michel ITMT1 Michel et al. ITMT1 
Liu et al. lTMfl Li et al. l2UU21 Li et al. lMMl Fox-Machado 



Table 2. Observed frequencies used in this work. The 
data were obtained from different STEPHI campaigns as 
detailed in the text. 



star v 

(MHz) 


star v 
( M Hz) 


HD 23628 191.8 
201.7 
376.6 


V624 Tau 242.9 
409.0 
413.5 
416.4 
451.7 
489.4 
529.1 


V650 Tau 197.2 
292.7 
333.1 
377.8 


HD 23194 533.6 
574.9 


V647 Tau 236.6 
304.7 
315.6 
374.4 
405.8 
444.1 


V534 Tau 234.2 
252.9 
307.6 
377.9 
379.0 
448.1 
525.0 



et al. I2002[) . The frequency peaks detected with a con- 
fidence level above 99% are summarized in Table El We 
note that the frequency resolution in a typical STEPHI 
campaign (three weeks) is Av ~ 0.5 /iHz. 

Figure ^ illustrates a colour-magnitude diagram of 
the Pleiades cluster in the region where the target stars 
are located. The filled circles corresponds to the S Scuti 
stars. The apparent magnitudes, V, colours, (B — V) and 
membership used are contained in the WEBDA database 
(Mermilliod 20013). The Pleiades shows differential red- 
dening with significant excess in the southwest with an 
average value of E{B — V) ss 0.04 (e.g. van Leeuwen 
IT<%)51 Breger ITMfl Hansen-Ruiz & van Leeuwen ITWH 
Pinsonneault et al. ^998). In particular, we adopted the 
reddening for individual stars as given by Breger (1986) 
if present, otherwise the average of E(B — V) = 0.04 was 
used. Stars known to be spectroscopic binaries were re- 
jected. 

Given the young age of the Pleiades cluster, we expect 
all the (5-Scuti stars to be at an early evolutionary stage on 
the main sequence. From the figure it follows that while 
V647 Tau and V624 Tau have similar masses, V650 Tau, 
HD 23194 and HD 23628 are slightly more massive. V534 
Tau, located near the red edge of the instability strip, is 
the coolest star in the ensemble. 



3. The modelling 

In this section we shall explain how we have computed 
the stellar models, the corresponding eigenfrequencies and 
how we have constrained the stellar parameters used in the 
forward comparison with the observed frequencies. 
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Table 1. Observational properties of the target stars 



Star 


HD 


ST 


niy 


B - V 


E(B - V) 


V sin i 


(3 














(kms -1 ) 






23628 


A7V 


7.66 ±0.03 


±0.211 ±0.005 


0.063 ± 0.008 


150 ± 15 


2.884 


V650 Tau 


23643 


A3V 


7.77 ±0.02 


±0.157 ±0.006 


0.027 ± 0.008 


175 ± 18 


2.862 




23194 


A5V 


8.06 ±0.02 


±0.202 ± 0.010 


0.076 ± 0.008 


20 ± 2 


2.881 


V624 Tau 


23156 


A7V 


8.23 ±0.01 


±0.250 ± 0.005 


0.046 ± 0.008 


20 ± 2 


2.823 


V647 Tau 


23607 


A7V 


8.26 ±0.01 


±0.255 ± 0.001 


0.057 ± 0.008 


20 ± 2 


2.841 


V534 Tau 


23567 


A9V 


8.28 ± 0.02 


±0.360 ± 0.001 


0.084 ± 0.008 


90 ±9 


2.788 



F 




•1 t 1 

0.0 0.1 0.2 0.3 0.4 

(B-V) 

Fig. 1. HR diagram of the Pleiades cluster. Only the re- 
gion around the instability strip is shown. The target stars 
are represented with filled circles. The blue and red edges 
of the instability strip are indicated by continuous lines 
and were taken from Rodriguez et al. (2000). An isochrone 
of 70 Myr is also shown. 

3.1. The stellar models 

We have considered stellar models with input physics ap- 
propriate to the mass range covered by S Scuti stars. In 
particular, the stellar models were computed using the 
CESAM evolutionary code (Morel ITW7| . The nuclear re- 
action rates are from Caughlan & Fowler l|1988[) . the equa- 
tion of state is from Eggleton et al. I|1973[) . the opaci- 
ties from the OPAL project (Iglesias & Roger 1996.) com- 
plemented at low temperatures by Kurucz data (Q991), 
and the atmosphere is computed using the Hopf's T(t) 
law (Mihalas 1978). The convection is described accord- 
ing to the classical mixing- length theory (MLT hereafter). 
In particular, we have considered an MLT parameter of 
omlt = I /Hp — 1-52, where H p is the pressure scale- 
height. We have considered a fixed value for «mlt because 
this parameter is not expected to have any significant in- 
fluence on the global structure of the evolution models 
considered for our target stars. 

We have considered different sets of initial chemical 
compositions in order to cover the Pleiades [Fe/H] ob- 
served range (see Section 0) and different initial He abun- 
dances. The values used are given in Tabled For each ini- 
tial chemical composition we have computed models with 



Table 3. Metallicity of the models. [Fe/H] = log(Z/X)» - 
log(Z/ X) , with (Z/X)q = 0.0245 from Grevesse & 
Noels itTMSl 



X 


Y 


Z 


[Fe/H] 


Q 


yv 


0.735 


0.250 


0.015 


-0.0794 


0.0- 


-0.2 


0.685 


0.300 


0.015 


-0.0488 


0.0- 


-0.2 


0.700 


0.280 


0.020 


0.0668 


0.0- 


-0.2 


0.725 


0.250 


0.025 


0.1484 


0.0- 


-0.2 


0.675 


0.300 


0.025 


0.1794 


0.0- 


-0.2 



and without core overshooting, ct ov = 0.20 in the former 
case (Schaller et al. fiMty . 

Although the final comparison with observations is 
done considering rotating models, our procedure requires 
the computation of non-rotating models and the cor- 
responding isochrones. Hence, sequences of non-rotating 
models were calculated for masses between 0.8 Mq and 5.0 
Mq , from the Z AMS to the sub-giant branch and the cor- 
responding isochrones computed for each sets of parame- 
ters in Table |3| In the following analysis we shall consider 
three ages, A, of 70 Myr, 100 Myr and 130 Myr. Finally, for 
comparison with the observations three distance moduli, 
d, of 5.50, 5.60 and 5.70 mag were considered except for 
models with [Fe/H]=— 0.0488, in which case an additional 
value of d = 5.39 mag was used. 

Figure ^ shows an example of such isochrones com- 
puted with the following parameters: [Fe/H]=0.0668, 
a ov = 0.2, A = 70 Myr and d — 5.70 mag. The isochrones 
were calibrated from [T e s, log(i/L Q )] to (B — V, Mv) by 
using the Schmidt-Kaler {19321 calibration for magni- 
tudes and the relationship between T e g and B — V of 
Sekiguchi & Fukugita (200TJ) for the colours. As can be 
seen in Fig. in the case illustrated here the isochrone 
matches the observed colour-magnitude diagram. For 
some combinations of high metallicity and low distance 
modulus the isochrone fit is not satisfactory but we have 
not rejected those combinations to allow for an indepen- 
dent determination of the cluster parameters from the os- 
cillation frequencies. Hence we have considered the entire 
combination of cluster parameters given in Table |21 and 
the three ages given above. 
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3.2. Photometric effect of rotation on the HR diagram 

It is known that fast rotation affects the position of a star 
in the colour-magnitude diagram (e.g. Tassoul 1978). In 
particular a ZAMS rotating model has a smaller luminos- 
ity and mean T g than a non-rotating model with the same 
mass and chemical composition. Also, the magnitude and 
the colour of a rotating star depend on the aspect angle, 
i, between the line of sight and the rotation axis. In par- 
ticular, a rotating star seen pole-on appears brighter and 
hotter than the same star seen equator-on. 

We now search for the range of masses and angular 
rotations suitable for each star. To do this we shall use 
the results of Perez Hernandez et al. (1999), in which the 
correction to the photometric magnitudes of non-rotating 
models needed to obtain those of rotating models were 
computed for main sequence stars of spectral type AO to 
F5. In this calculation, we shall take into account the ob- 
served v sin? for each star. As an example, Fig.[2illustrates 
the corrections due to rotation in the particular case of the 
S Scuti type star V650 Tau with wsini = 175 km s _1 . The 
actual position of the star is shown with a thick dot upon 
which there is a cross that gives the errors in the esti- 
mation of (Mv)o and (B — V)o calculated according the 
following expressions: 

*[(B - V) ] = V° 2 (B -V)+ <r 2 [E(B - V)}, (1) 
a[(M v )o] = V° 2 (d) + <? 2 (rn v ) + (3.2a[E(B - V)]f , (2) 

where <j(d) = 0.06 is the error in the distance modulus. 
The errors in the magnitude My, colour index (B — V) 
and reddening E(B — V) are given in Tabled 

The range of masses and rotational velocities of the 
star depends on the stellar parameters considered and 
hence needs to be computed for each combination in 
Tabled In particular, in the example illustrated in Fig. [21 
the same cluster parameters as in Fig. ^ were considered 
([Fe/H]=0.0668, a ov = 0.2, A = 70 Myr and d = 5.70 
mag). The isochrone is represented by a thick continu- 
ous line and some of its models are shown by dots. For 
each model on the isochrone there is a small track corre- 
sponding to the photometric corrections due to rotation 
computed by using the projected velocity usini of the 
star and changing the angle of inclination from i = 90° to 
the minimum angle i = £ m in corresponding to ~ 0.90f2b r , 
where fibr is the break-up angular rotational velocity. This 
angular velocity is given by the balance between the grav- 
itational attraction and the centrifugal force at the equa- 
tor. The dotted lines give the corrections at fixed angles: 
i = 90° (the closest line to the isochrone), i = 50° (inter- 
mediate line) and i = 39° (the reddest line). It is evident 
that the actual position of the star, including the error 
box, is associated to a domain of non-rotating models. In 
this case we have obtained for the star considered a range 
of masses of [1.86, 1.91]Af Q and a range of angles of in- 
clination of [39°, 50°]. Nonetheless, in order to take into 
account possible errors coming from the isochrone cali- 
bration we shall consider a wider range of aspect angles 



M/M a =2.45 




3 5 I i I i i i I i i i I i I 

0.2 0.4 

B-V 

Fig. 2. Colour-magnitude diagram illustrating the photo- 
metric corrections due to rotation in the case of V650 Tau 
(vsini = 175km s~ x ). The thick continuous line is the 
same isochrone represented in Fig. ^ ([Fe/H] = 0.0668, 
a ov = 0.2, A = 70 Myr, d — 5.7 mag). Some models are 
shown with filled circles upon the isochrone. The masses 
of the first and last models are indicated. The small tracks 
represent the photometric corrections due to rotation for 
each the models indicated on the isochrone when i runs 
from i = 90° to % = i m i n - The corrections at fixed i are 
represented by dotted lines. 



i = [imin) 90°]. The same procedure is carried out for the 
other stars. Also one needs to consider all the metallici- 
ties, ages, overshooting parameters and distance moduli 
for the whole set of stars. 

Having obtained ranges of M and values of Win, we 
proceed to find the corresponding range of angular rota- 
tional velocities f2, for which an estimate of the equato- 
rial radius, i? C q, of the rotating model is needed. Under 
the assumption of uniform rotation and approximating 
the surface of the star by a Roche surface (see e.g. Perez 
Hernandez et al. I1999f) one has: 

3 f 1 r nl 

R cq ~ R Q — cos < - [tt + arcoswj > , (3) 

where ui ~ fi/fi br , f7g q ~ 8GM/(27i?g) and R is the 
polar radius. As noted in Perez Hernandez et al. (^999), 
the polar radius can be approximated by the radius of 
a spherical non-rotating model with the same mass and 
evolutionary state (since our stars are close to the ZAMS, 
a non-rotating model of the same mass and age is suitable 
for this error-box estimation). 
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1.90 
1.85 
1.80 
1.75 
1.70 
1.65 
1.60 



0.0 0.5 1.0 1.5 

(10~ 4 rad/s) 



2.0 



Fig. 3. Range of rotation rates against masses estimated 
for the six stars from the photometric corrections due 
to rotation applied to the same isochrone depicted in 
Figure [3 The vertical error bars give the estimated range 
of mass for each star. The crosses give the values of the 
break-up angular rotational velocity, f^br- The ranges of 
ft and M are associated with stars as follows: dot-dashed 
line for V650 Tau (wsini = 175 km s _1 ); long-dashed line 
for HD 23628 (wsini = 150km s" 1 ); three-dots dashed 
line for HD 23194 (wsini — 20 km s _1 ); dashed-line for 
V647 Tau (vsini = 20 km s _1 ); continuous line for V624 
Tau (vsini = 20km s _1 ) and dotted line for V534Tau 
(fsini = 90 km s _1 ). 



Since on the other hand the angular rotation is related 
to the equatorial radius by 



n 



(v sini) 



bs 



R 



(4) 



cq 1 



it is possible to carry out a simple iterative process to 
obtain a range of possible rotations, [O m ; n , fi max ], from 
the previously obtained range of angle of inclinations 
[*min,90°] and masses. Since for given data parameters, 
the mass of the star is obtained with an uncertainty of 
~ 0.02 Mq this iterative process was carried out for just 
one mass. 

As an example, Figure |31 shows the ranges of fl ob- 
tained for the six stars as a function of the mass of 
the models considering the same cluster parameters as in 
Figured Lines of different types give the estimates for the 
different stars as indicated in the figure. A similar proce- 
dure needs to be considered for all sets of parameters in 
Table El 

With the information obtained above we compute evo- 
lutionary sequences of rotating models with the same in- 
put physics as the non-rotating models but with appropri- 
ate initial angular rotational velocities in order to match 
as closely as possible the extreme values in the interval 
[^min, ^max] at the final age. Solid-body rotation and con- 
servation of global angular momentum during the evolu- 
tion were assumed during the calculus. In order to have 



a better overview of the rotation rates of the stars we 
have also computed models with an intermediate value of 
^mid ~ 0.5fibr- A total amount of 1620 sequences of rotat- 
ing models were finally obtained for the whole ensemble 
of stars. 



3.3. Theoretical oscillation frequencies 

The eigenfrequencies of the rotating models previously de- 
scribed have been calculated using the oscillation code 
FILOU (Tran Minh et al. HHEI Suarez \W)2] . We have 
considered frequency perturbations up to second order 
in the rotation rate. Frequencies of the modes are la- 
belled in the usual way: n, I, m for the radial order, degree 
and azimuthal order. The radial orders of a given mode 
are assigned according to the Scuflaire criterion (1974) 
with n > for the p modes, n < for the g modes, 
n = 1 for the fundamental radial mode of degree and 
1, and n = for the fundamental mode with the 1 = 2. 
Eigenfrequencies were computed up to I = 2. For geomet- 
rical reasons higher degree modes are expected to have 
considerably smaller amplitudes. The theoretical frequen- 
cies cover the frequency range of the observed pulsation 
peaks (see Table |2J)- Coupling between modes is not con- 
sidered in the present work. 

The estimated interval of fi ro t for all the stars may be 
as large as Ar2 rot ~ 1.6 x 10 -4 rad/s which in terms of 
cyclic rotational frequency corresponds to Av m t ~ 25 /iHz 
(see for instance Fig. 0J. After some tests, we found that 
a satisfactory comparison between the observed and the- 
oretical frequencies cannot be achieved by using only the 
frequencies computed for models with three f2 ro t within 
this interval. To overcome this difficulty once we had the 
mode frequencies for each series of models (three f2 ro t for 
fixed M, [Fe/H], «mlt, oe ov , d and ^4), we proceeded to 
interpolate the results on fl mt in 21 equally spaced points 
covering the interval [T2 m in, !J m ax]- A quadratic spline in- 
terpolation was applied. In order to evaluate the goodness 
of the interpolation we also computed the eigenfrequencies 
for several randomly selected models. We have found that 
there is a good agreement between the interpolated and 
theoretical frequencies up to v ~ 600 /iHz (n ~ 6), the av- 
erage of absolute separation being \v\ lx \_ — v ca \ | - 0.3 /iHz. 
The disagreement found beyond this value can be ex- 
plained by the fact that the second order effect of rotation 
is enhanced at higher frequencies. In any case, as can be 
seen in Tabled the highest frequency we are dealing with 
is ~ 574 (Mz. 



4. Method of comparison 

Once we have the mode frequencies for all sets of parame- 
ters we compare the observational frequencies {v hs) with 
the theoretical ones (f ca i) at each interpolated f2 ro t. 

Let us consider a rotating model with given parameters 
for just one star. Then, for all the possible combinations 
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between the observed and computed frequencies, we com- 
pute the quantity x 2 given by 



1 N 



(5) 



where N is the number of observational frequencies. In 
this computation it is understood that each theoretical 
frequency is assigned to one observed frequency at most. 

As a first step, in order fully to exploit the collective 
behaviour of stars within open cluster we consider the so- 
lutions of the six stars with given cluster parameters of 
[Fe/H], a ov , d and A. The models corresponding to each 
star differ in the mass and the rotation rate. As an ex- 
ample, Fig^ shows x 2 against angular rotational velocity, 
tt, for the six models corresponding to the same parame- 
ters as in Fig.Elbut without overshooting ([Fe/H]=0.0668, 
a ov = 0.0, g?=5.70 mag and A=70 Myr). Each panel corre- 
sponds to the star indicated in the figure. For clarity, only 
solutions with x 2 < 20 (/iHz 2 ) are shown. Similar features 
are found for other metallicity, distance moduli and ages, 
but the lowest \ 2 may appear with rather high values. 
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0.0 



0.1 



0.2 



[Fe/H] 



Fig. 5. Minimum values of the mean square root of the 
difference between observed and theoretical frequencies 
for models with the same values of [Fe/H], a OVl d and 
A against [Fe/H]. For clarity, only models without over- 
shooting are shown. The symbols are related to different 
values of A and d as follows: asterisks (d — 5.39 mag), 
squares (d = 5.50 mag), diamonds (d = 5.60 mag), circles 
(d = 5.70 mag), small (70 Myr), middle (100 Myr) and 
big (130 Myr). 




0.6 0.8 1.0 1.2 
(x10~ 4 rad/s) 



Fig. 4. x 2 against Q for six models with common param- 
eters [Fe/H], a ov , d and A. 



Since we expect the same [Fe/H], d and A for all the 
stars in the cluster we shall assume that the best fits 
should happen simultaneously in the six stars despite dif- 
ferences in rotation rate and masses of the models. Hence 
we proceed to calculate the mean square root of the lowest 
X 2 found in each model for given common parameters by 
means of the following expression: 



\ 



1 6 

i— 1 



(6) 



Figure shows e against [Fe/H] for models without over- 
shooting. For the sake of clarity only the points with 
e < 3.0 /iHz are shown. It can be seen that the solution 
with the smallest e corresponds to [Fe/H]=0.0668. A plot 
of similar characteristics is obtained for models with over- 
shooting. 

Since only those identifications obtained from eq. (5) 
with low x 2 values are of interest, we require that for each 
star the solution must have x ^ 3.5 /iHz. In Fig. we 
reject those solutions that have at least one model with 
X > 3.5 /iHz for all the parameters and combinations be- 
tween observed and theoretical frequencies. With this re- 
striction we found that only those solutions associated 
with the models represented by the four lowest points in 
Fig. should be considered. Similar solutions were found 
for models with overshooting. All these models are listed 
in Tableland will be analysed in detail later. 

We shall use a geometrical argument to reduce fur- 
ther the number of solutions so far available for each star. 
To this end, we take into account that the visibility of 



L. Fox Machado et al.: A seismological analysis of S Scuti stars in the Pleiades cluster 



7 



Table 4. Ensembles of models with and without over- 
shooting, which, after applying a threshold of x < 3.5 /iHz, 
have potential valid solution (see text for details). 







[Fe/H] 


= 0.0668 






Qov = 


0.00 




Otov = 


0.20 




d 


A 




d 


A 




(mag 


(Myr) 




(mag 


) (Myr) 


CI 


5.70 


70 


Dl 


5.70 


70 


C2 


5.70 


130 


D2 


5.60 


100 


C3 


5.60 


70 


D3 


5.60 


130 


C4 


5.60 


100 


D4 


5.70 


100 



0.6 



.. i-o 






*•••., 1=1, m=0 


1=1, |m[=f - 




1=2, m=0 










t=_2 L | mj = 1 








' ' 1=2, | m| =g S N - ' 

x x 






- r- - r " 1 





20 40 60 80 

i (degree) 

Fig. 6. Variation of visibility amplitude with the inclina- 
tion angle for radial I = (continuous line) and non-radial 
I = 1 (dotted lines) and I = 2 (dashed lines) modes. 



each mode depends on the angle of inclination. Following 
Pesnell (1985), we illustrate in Fig.^the spatial response 
function Si m against i, for mode degree I = 0,1,2. For 
simplicity, limb-darkening has been neglected. It can be 
seen that for i « 90° only even I + m values can be de- 
tected, while for i « 0° only modes with m = will be 
visible. 

We introduce an hypothesis based on the visibility of 
the modes: we reject any solution with one or more modes 
with S[ < 0.1. After applying this constraint we were 
able to limit the analysis to four sets in Table 0] two with 
a m = 0.0: CI (d = 5.70, A = 70), C3 (d = 5.60, A = 70) 
and two with a ov = 0.2: Dl (d = 5.70, A = 70), D2 
(d = 5.60, A = 100). We note that while for the stars 
HD 23628 and HD 23194 only a few identifications (< 8) 
remain, the number of identifications for V624 Tau, V534 
Tau, V647 Tau and V650 Tau remained larger than 100 
in most cases. 



5. Results and discussion 

In order to discuss the results we shall introduce a param- 
eter A associated with each star in each identification and 
given by 

A = max( | v ohs - t'cai | ) • (7) 

Table|S]lists the number of solutions for the six stars in 
each ensemble. Different maximum values of A have been 
considered. Certain features can be derived from these 
general results: 

— The resulting masses for each star are similar for all 
the solutions. 

— The number of solutions for a given star at fixed A 
associated with models with overshooting is smaller 
than without it. 

— There are no solutions with A < 3.0 /iHz for the stars 
V534 Tau and V647 Tau, while the star HD 23194 
shows solutions even with A < 1.0/iHz. 

— In order to find a set of models that has at least one 
solution simultaneously for all the stars a value of A ~ 
3.5 /iHz is needed. This solution is found for CI that 
corresponds to the cluster parameters [Fe/H] =0.0668, 
a ov = 0.0, d = 5.70 and A = 70, and it coincides with 
the minimum in Figure [3] 

— The remaining ensembles in Table |S] have at least 
one solution simultaneously when the value of A is 
increased slightly. In particular, for C3 a value of 
A ~ 4.5 /iHz is needed to obtain a solution while for 
Dl and D2 a value of A ~ 3.8 /iHz is required. 

We have found that the identifications and the range 
of aspect angles derived for each star are similar for all the 
solutions. This agreement may be due to the fact that the 
estimated range of mass and radius, and hence of mean 
densities, is similar for all the solutions. At the ages consid- 
ered the overshooting has a negligible effect on the models. 

Table El summarizes the parameters estimated for the 
cluster as well as for each star associated with the possible 
identifications. The radial orders of these modes are con- 
sistent with growth rate predictions (Suarez et al. [2002 ) . 
While these identifications cannot be considered as defini- 
tive but rather as "best fit solutions" , some information 
on the stars could be derived. On the one hand, those stars 
with smaller masses (V624 Tau, V534 Tau and V647 Tau) 
have more than six frequencies and could have simultane- 
ously excited radial and non-radial oscillations although at 
most two radial modes. Following Hernandez et al. fl998), 
at most two radial modes were found to be excited for 
the 5 Scuti stars in the Praesepe cluster. On the other 
hand, the most massive stars, V650 Tau, HD 23628 and 
HD 23194, with a mean mass of 1.83 Mq present on aver- 
age few observational frequencies (N < 4) that might be 
associated only with non-radial modes. With the present 
observational data it remains unclear, however, whether 
this is a general trend for <5 Scuti stars. 
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Table 5. Number of possible solutions for models in each ensemble for different values of A (in /xHz). 
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V624 Tau 


1.70 








12 
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1.86 





2 
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1.68 
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1.69 
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1.69 
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1.86 














1.86 











3 


HD 23194 
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1 
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1.79 


1 


1 
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HD 23668 


1.82 








3 


4 


1.84 





2 


3 
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A comparison of observed and computed frequencies 
is shown in Figure In each panel the solid lines repre- 
sent the observational frequencies with their correspond- 
ing amplitudes. The "best fit" theoretical frequencies are 
represented by dotted lines. The stars are arranged from 
top to bottom according to their magnitude. 

From this figure it follows that there is an asymmetric 
triplet in the observational spectrum of V624 Tau. We 
have found that for all the solutions the identification is 
always a (I = 1, n = 3) mode split by rotation. In turn, 
this implies that an independent estimate of the mean 
rotation rate can be obtained from the rotational splitting, 
up to second order, given by 

We have found complete agreement between the cyclic ro- 
tational frequency computed from the above expression 
and those derived from the mode identification, the differ- 
ences being ~ 0.03 — 0.04 //Hz. We have also found a fairly 
good agreement between the observations and the models 
in the differences v- m — vq and v +m — vq. 

If the estimated inclination angles for the stars are 
correct, all of them could be rapid rotators (v mt > 10/xHz) 
except for V624 Tau with v rot < 5^Hz. HD 23194 and 
V647 Tau, whose projected rotational velocities are low, 
v sin i — 20 kms -1 , would have equatorial velocities as 
high as 90 and 70 kms -1 respectively. 

In the present work the distance is found to be my — 
My = 5.60 — 5.70 mag (132 — 138 pc) and agrees very well 
with that derived recently by independent determinations 
using independent techniques and data (Pan et al. |J004 



Munari et al. I2tjfjl Soderblom et al. l2"tjfJ5l Percival et al. 
2005). However, although the isochrones computed with 
sub-solar metallicity (Z = 0.015 and Y = 0.30) match well 
the observational data with a distance modulus of d = 5.39 
mag in agreement with that of HIPPARCOS within 1-cr 
error, the resulting solutions (shown by asterisks in Fig.[5j 
do not lead to good fits. Moreover, the plot of e against d 
which we do not reproduce here, shows better fits for fur- 
ther distances. For solar metallicity, in order to reproduce 
the HIPPARCOS MS of the Pleiades, unrealistically high 
helium abundances are required (e.g. Y — 0.34 Belikov et 
al.UHS Y = 0.37 Pinsonneault et al. H555^ . 

6. Conclusions 

In this study we have performed a seismological analysis of 
six 5 Scuti stars belonging to the Pleiades cluster to iden- 
tify their frequency modes. To the best of our knowledge 
this group of variables constitutes the most statistically 
significant sample of S Scuti stars analysed in an open 
cluster to date. 

Rotational effects were considered in different stages of 
the analysis: first when determining the star positions in 
the HR diagram and second when computing stellar mod- 
els that approximately reproduce the evolutionary stage 
of the stars, and finally when computing theoretical oscil- 
lation frequencies in order to construct seismic models for 
target stars. The comparison between observational and 
computed frequencies was carried out by a least-squares 
fit. There is a large number of possible solutions partly 
because neither the equatorial velocity nor the inclination 
angle are known a priori. In order to limit the number 
of possible solutions we used the relationship between the 
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Table 6. Summary of parameters derived for the Pleiades cluster as well as for S Scuti stars with the possible 
identification of the frequency modes. 



[Fe/H] = 0.0668, aw 


= [0.00-0.20], mv - 


M v =[5.60-5.70], A =\ 


70-100] x 10 
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Identification 
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V 
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(uHz) 


(n, I, m) 




( M Hz) 


(n, i, m) 


V624 Tau 


242.9 


(1,0,0) 


V650 Tau 


197.2 


(1,1,1) 


M = [1.68-1.72] M 


409.0 


(3,1,1) 


M = [1.84-1.88] M Q 


292.7 


(0,2,-2), (2,2,2) 


frot = [3-5] fiRz 


413.5 


(3,1,0) 


u rot = [25-28] ^Hz 


333.1 


(3,1,1), (3,2,2) 


i = [37°-67°] 


416.4 


(3,1,-1) 


i = [60°-70°] 


377.8 


(2,2,-2), (3,1,0) 




451.7 


(3,2,-2), (4,0,0) 






(3,1,-1), (3,2,1) 




489.4 


(4,1,0),(4,1,1) 
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(4,2,-1), (4,2,-2) 












(5,0,0) 








V534 Tau 


234.2 


(1,1,1) 

V 3 ? / 


HD 23628 


191.8 


(0,2,2) 


M = [1.65-1.69] M© 


252.9 


(1,1,0) 


M = [1.82-1.86] M 


201.7 


(1,1,1) 


i/rot = [14-16] /iHz 


307.6 


(2,0,0), (2,1,1) 


u Iot = [24-26] ^Hz 


376.6 


(2,2,-2) 


i = [59°-79°] 


377.9 


(2,2-l),(3,0,0) 


i = [53°-59°] 








379.0 


(2,2-l),(3,0,0) 












(3,1,1) 










448.1 


(3,2,-l),(4,0,0) 










525.0 










V647 Tau 


236.6 


(1.1,1) 


HD 23194 


533.6 


(5,1,1), (5,1,-1) 


M = [1.68-1.72] M 


304.7 


(1,2,1) 


M = [1.78-1.82] M Q 


574.9 


(5,2,0),(5,2,1) 


frot =10—11 /iHz 


315.6 


(2,0,0), (2,1,1) 


v I0 t = [14-23] /xHz 






i = 17°-18° 


374.4 


(2,2,-1) 


i = 7°-12° 








405.8 


(3,1,0) 










444.1 


(3,2,0) 









amplitude visibility, Si m and the aspect angle, i. As a re- 
sult we found few solutions for each star, suggesting the 
existence of only p modes of low degree, low radial order in 
all the stars. For the less massive stars, solutions with at 
most two radial modes were also possible. These solutions 
have associated ranges of stellar parameters for each star. 
Most of the stars could be rapid rotators according to the 
estimated angle of inclination, i. 

The best fits between observational and theoretical 
frequencies are achieved for global cluster parameters of 
[Fe/H] = 0.0668 (Z = 0.02, X = 0.70), A = 70-100Myr 
and d = 5.60 — 5.70 mag. The derived distance modulus for 
the Pleiades cluster agrees with that of the main sequence 
fitting method, in spite of the fact that the isochrones with 
sub-solar metallicity closely matches the Pleiades HR di- 
agram with the HIPPARCOS distance. 
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